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We examine the contributions from triple excitation cluster operators in the relativistic coupled- 
cluster theory in atoms and ions. For this, we propose a tensor representation of the triple cluster 
operator. Based on this representaion and using diagrammatic analysis, we derive the linearized 
coupled-cluster equations for single, double and triple excitation cluster operators. The contributions 
from the triple cluster operators to the hyper fine structure constants of single-valence systems Rb 
and Sr + are analysed using the perturbed triples. 



I. INTRODUCTION 

The coupled-cluster (CC) theory [IH3] is an all-order 
many-body theory. It has proved to be one of the most 
reliable and accurate methods for precision atomic the- 
ory calculations. Apart from atomic systems [H[5], it has 
also been used with great success in nuclear [5] , molecular 
[7] and condensed matter [5] calculations. In atoms and 
ions, calculations using the relativistic coupled-cluster 
(RCC) theory has provided some of the best results. 
These include calculation of atomic electric dipole mo- 
ments |H [5] , hyperfine structure constants [25] , elec- 
tromagnetic transition properties [HI 02] and most im- 
portantly the NSI-PNC [T3]. 

Among different flavours of CC theory, the coupled- 
cluster singles and doubles (CCSD) is a widely used ap- 
proximation. However, for precision atomic calculations 
it is imperative to estimate the contributions from clus- 
ters of higher excitations. For CCSD approximation, the 
triple excitation cluster operators is the closest level of 
excitation neglected in the calculations. So, the leading 
order correction to CCSD is the effect of triple excitation 
cluster operators. Further more, the triple excitations 
are expected to have significant contributions in open 
shell systems. Due to the N„ scaling, where N v is the 
number of virtual orbitals, inclusion of triple excitation 
cluster operators pose severe computational challenges. 
A practical approach is selective inclusion of triple exci- 
tation cluster operators. Such calculations are important 
to make uncertainty estimates. In this work we examine 
contributions from the valence triple excitation cluster 
operators S3 in RCC and propose a representation of S3. 
Furthermore, to quantify the importance we carry out ex- 
tensive calculations with different forms of perturbative 
triple excitations. 

Atomic parity non-conservation (PNC) is one class of 
atomic theory calculations, where precision theory cal- 
culations are important and uncertainty estimates are a 
must. The atomic theory results of PNC observable when 
combined with the experimental data provide estimates 
of parameters in standard model (SM) of particle physics 
[13] . Any deviation from the predictions of SM is an indi- 
cation of new physics. The PNC in atoms occurs in two 
forms, nuclear spin-independent (NSI) and nuclear spin- 
dependent (NSD). The former, has been theoretically and 



experimentally studied in great detail, and the most ac- 
curate theoretical [14] and experimental [15] results are 
in the case atomic Cs. For the later, however, there are 
few theoretical studies. These are using many-body per- 
turbation theory (MBPT) [TB], configuration interaction 
(CI) [mm] and CI+MBPT [T9J [20]. 

We recently proposed an RCC based method to incor- 
porate nuclear spin-dependent interaction Hamiltonian 
as perturbation. The method is used to calculate the 
NSD-PNC of Cs, Ba+ and Ra+ [21] with associated rms 
uncertainties of 7%, 4.4% and 7.6%, respectively. The 
details of the proposed theory are presented in another 
work of ours [22]. We believe that its possible to reduce 
the uncertainty, and the first step towards this could be 
the inclusion of triples cluster operators in RCC. 

The paper is divided into nine sections. In Section. 
II, we give a brief review of CCSD. It is based on our 
previous works (23] [24] on RCC of closed-shell and one- 
valence systems. Then the next section, Section. Ill, 
forms the core of the present work and describes the per- 
turbative S3. It discusses the possible chanels through 
which S3 can arise and describes the tensor structure. 
In Section. IV, we give linearized RCC equations for 
singles, doubels and triples in terms of the CC excita- 
tion amplitudes. The HFS constants calculation using 
CCSD is breifly demonstrated in Section. V for the easy 
reference. A detailed description of HFS terms in RCC 
properties calculations and diagrams from the perturba- 
tive triples are given in the Sections. VI and VII. And in 
Section. VIII, we present and discuss our results. 

II. BRIEF REVIEW OF RCC IN CCSD 
APPROXIMATION 

The Dirac-Coulomb Hamiltonian which accounts for 
the leading order relativistic effects of atoms or ions with 
N electrons is 

N 1 

h dc = J2 [can ■ p» + (A - i)c 2 - V N {n)] + J2 — , 

i— 1 i<j 3 

where oli and /3 are the Dirac matrices, p is the linear 
momentum, Vjv(r) is the nuclear Coulomb potential and 
last term is the electron-electron Coulomb interactions. 



For one-valence systems it satisfies the eigen value equa- 
tion 



(2) 



where and E v are the atomic state and energy, re- 
spectively. In the CC method, the is expressed in 
terms of T and S, the closed-shell and valence cluster 
operators respectively, as 



(3) 



where is the one- valence Dirac-Fock reference state. 
It is obtained by adding an electron to the closed-shell 
reference state, \$> v ) = aJ,|$o). For an N electron system, 
which may be atom or ion, the cluster operators are 



N-l 



1=1 

N 



(4a) 
(4b) 



The index i represents the level of excitation (loe) of the 
cluster operators. Note that for T loe is allowed up to 
the N — 1 core electrons, where as for S it is up to N as 
it includes the valence electron. One major impediment 
to a full scale CC calculation is, the number of cluster 
operators proliferates exponentially with i and calcula- 
tions are unmanageable beyond the first few loe. This 
difficulty, as such, does not diminish the applicability of 
CC. Most dominant correlation effects are incorporated 
in the first few loe and the approximation referred to as 
the CC singles and doubles (CCSD) provides a very good 
description of the many-body effects. In this approxima- 
tion 



T = T 1 +T 2 , and S = Si + S». 
The operators in second quantization notation are 

Ti = E *S4°»' and T2 = 4 E €l a l a l a b a a- 



Si = y^sga^q^, and S 2 



abpq 



^ sf, 9 „aiata a a v . 



apq 



va p q 



(5) 



(6) 



(7) 



Here, t".\ and s\\\ are the cluster amplitudes. The indexes 
abc... (pqr . . .) represent occupied (virtual) states and 
vwx . . . represent valence states. The operators T\ {S\ 
) and T 2 (S 2 ) give single and double replacements after 
operating on the closed-(open-)shell reference states. The 
diagrammatic representations of T and S are shown in 
the Fig. [1] 

The closed-shell CC operators are the solutions of the 
nonlinear coupled equations 



«|i? N |$0> =0, 



(8a) 
(8b) 
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T 3 
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FIG. 1. Diagrammatic representation of Ti and T2 operators. 
The orbital lines with up (down) arrow indicate particle (hole) 
states. 



where /Tn = e T H^e T is the similarity transformed 



Hamiltonian and the normal order Hamiltonian = 
H-{$o\H\$ ). The states and |$^) are the singly 
and doubly excited determinants, respectively. The de- 
tails of the derivation are given in Ref. [23 . The one- 
valence CC operators on the hand are obtained from the 
solutions of the equations 

($p\H n + {H^S}\$ v ) = £5 tt <$5|Si|$„), (9a) 

(9%\H N + {h7s}\$ v ) = E^iS&WQv), (9b) 

where = E v — Eq, is the attachment energy of the 
valence electron. The details of the derivation of Eq. ^ 
we provide in our previous work |24) . 



III. PERTURB ATIVE TRIPLES IN RCC 

Inclusion of perturbative triples to the CCSD approx- 
imation is refereed to as the CCSD(T) approximation. 
Within this approximation, Eq. (4b) is 



T = T 1 +T 2 +T 3 and S = Si + S* 2 + S 3 . 



(10) 



Here, T3 and S3 are the perturbative core and valence 
triple excitation cluster operators, respectively. Like the 
single and double operators, second quantized form of 
the triple excitation cluster operators are 



Ts = 31 E tlZ a t a l a t a c a l> a a> 



abcpqr 



Ss = V\ E S l q ab a i a \ a i a b a aa v - 



(11a) 



(lib) 



abpqr 



Previous calculations have shown contribution from T 3 to 
the properties of one-valence systems are much smaller 
than S3, which is evident from the previous calculations. 
In particular the results reported in ref. |25j , where RCC 
is used. For this reason, in the present work, we consider 
and analyze in detail the contribution from the 5*3 cluster 
operators. The S3 operator in RCC arise from two chan- 
nels of contractions. First, residual Coulomb interaction 
(V 2 ) perturbs the open shell operator S 2 , and second, 
V 2 perturbs the closed-shell operator T 2 . Although, the 



contributions from the triples to properties are small, it 
is imperative to include for high precision calculations. 
The calculations related to discrete symmetry violations 
in atoms and ions belong to the class which require high 
precision. The need is even higher for open shell systems. 




(a) 




(b) 



(c) 



A. Triples from S2 



FIG. 3. Perturbative triple cluster operators (S3) diagrams 
arising from the term V2T2. The dashed line here is to repre- 
sent the Coulomb interaction. 



The Goldstone many-body diagrams of the first chan- 
nel, 52 perturbed triples, are shown in Fig. [2j There are 
three topologically distinct diagrams. The first two arise 
from the contraction of V2 with a virtual orbital of S2 
and the contribution is 



j>qr _ 

vab(vs) 



1 

pqr 



1 r „ps 
sb va 



J sa°vbJ 1 



(12) 



where, the energy denominator e 



pqr 
vab 



Cf,, with e as the orbital energies and the matrix 



elements in general are v. 



(kl\l/r 12 \ij), and 



(kl\S-z\ij). The subscripts (vs) indicate the cluster is from 
the contraction with S2 through a virtual orbital. In a 
similar way, the last diagram in Fig. [2] arises from the 




(a) 



(b) 



(0) 



FIG. 2. Perturbative triple cluster operators (S3) diagrams 

arising from the term V2S2. The dashed line represents the 
Coulomb interaction. 



contraction of the core orbital and the contribution is 

-1 



vab(cs) 



E v ar pq 
u cb J vc 

^vab r 



(13) 



B. Triples from T 2 

Like in the S2 perturbed triples, there are three Gold- 
stone diagrams which contribute to the T 2 perturbed 
triples and these are shown in Fig. [3] In the figure, 
the last diagram arises from the contraction of a virtual 
orbital and contribution is 



s pqr 

vab(vt) 



pqr / j u vs u abi 
vab r 



(15) 



where, the subscript (vt) indicates perturbed T 2 and con- 
traction of virtual orbital. The other two diagrams in the 
figure arise from the contraction of core orbital and con- 
tribution is 



vab(ct) 



= — V (v pc t 

pqr / j \ u va u 



qr + v cq t pr ) 

cb ' u va i cb> 



(16) 



The subscript (ct) indicates perturbed T 2 and contraction 
of core orbital. One key difference is noted when the 
above expressions are compared with the S2 perturbed 
triples. The number of diagrams arising from the core 
and virtual orbital contractions are interchanged in the 
two cases. The T 2 perturbed triples is then 



pqr pqr 

S vab(t) ~ S vab(vt) 



„pqr 

vab(ct) 



(17) 



The total perturbed triples is the sum of the two contri- 
butions 



The negative sign follows from the application of Wick's 
theorem in the operator contractions. It is also evident 
from the rules of Goldstone diagram evaluation [26] . Ac- 
cording to which the phase of a diagram is {-l) l+h , where 
I is the number of loops and h is the number of internal 
core lines. For the present case, the diagram in Fig. [2p 
has one internal core line and no loops. The subscript 
(cs), like in previous case, indicate the origin of the term. 
Collecting the terms, the S2 perturbed triples is 



vab 



vab(s 



't)ofe(t)' 



(18) 



All together, there are six perturbative S3 diagrams. 
Three each from the S2 and T 2 perturbations. In con- 
trast, for the perturbative T3 there are only two diagrams 
and one channel, T 2 perturbation. 



C. Tensor structure of S3 



s P1 r _ „PQ r 
vab(s) vab(vs) 



The two component 



and 



j>qr 

vab(cs) ' 



vab(c 



(14) 



have different 



number of terms as 5* 2 is topologically asymmetric. 
Closed-shell triples T 3 , on the other hand, have one term 
each. 



The diagrammatic form of triple excitation cluster op- 
erators as shown in Fig. [T]are convenient representations. 
However, it runs into serious difficulties while decompos- 
ing into angular and radial parts. The central vertex, 
consisting of four lines, has no viable equivalent tensor 
representation. To arrive at a consistent representation 



4 



of the tensor structure of the triple cluster operator, we 
analyze the angular reduction of perturbative triple exci- 
tation diagrams shown in Fig. [2] and [3] As an example, 
we examine the perturbative triples diagram in Fig. [2ja). 
The angular reduction of the diagram into phase factor, 





FIG. 4. Angular reduction of the perturbative S3 diagram 
shown in Fig. ^ a). 



6j-symbol and an irreducible free angular diagram are 
shown in Fig. |4j The later represents the combinations 
of 3j-symbols to represent geometric part of the matrix 
element and remaining is the physical part, this follows 
from the Wigner-Eckert theorem. 

There are other diagrammatic representations of the 
tensor structure of S3. An example is the one given in 
Ref. |27j , where an intermediate line is coupling of total 
angular momentum and tensor operator. In the present 
case, spin-orbitals are coupled pairwise to represent a 
matrix element of tensor operators of rank fcj (i = 1, 2, 3) 
which are again coupled. The tensor representation of S3 
in explicit form is 

S 3 = ]T s 3 (k 1 ,k 2 ,k 3 ){C kl (r 1 )C k2 (r 2 )} k3 

ki,k2,k3 

xC fe3 (f 3 ), (19) 

where C^r) are c-tensor operators of rank fc, and ?i are 
unit vectors in the coordinates of the ith electron. The 
notation {. . .} fca indicates coupling of two c-tensors to a 
k 3 ranked c-tensor. The angular momenta and the rank 
of the tensor operators must satisfy the triangular con- 
ditions \j v - j p \ < fel < j v + j p , \j a - j q \ fc 2 «S ja + jq, 

\jb~jr\ < h «S jb+jr, and |fci-fc 2 | < k 3 sC kx+k 2 . At the 
same time, spin-orbitals must satisfy the parity selection 
rule (-1) 1 »+ 1 p = (-iy«+h = (_i)'b+ir. One important 
property of the representation considered here is, the ver- 
tices in the tensor form of the triple operator can be inter 
changed with appropriate phase factor. In other words, 
there is an inherent symmetry in the coupling sequence 
considered. 



there is a large increase in the number of diagrams. An 
approximation, which incorporates the leading order ef- 
fects of triple cluster amplitude but with much less com- 
putational complexity is the linearized treatment of the 
triples. The number of cluster amplitudes, however, are 
still large. For this we consider the inclusion of valence 
triples S3 in the linearised RCC. From the definition of 
i?N, introduced in Eq. pi), linearised RCC is equivalent 
to the approximations 



(20a) 
(20b) 



To analyze the contributions from 6*3, the valence cluster 
operator S = Si + S 2 + S3, however, T = T 1 +T 2 . The 
RCC equations of the single and double excitation cluster 
amplitude, Eq. ([9]), in linear approximation are 

($S|£T N + {H^T} + {jSnS}|$„> 

^ =E: tt {^ v \S 1 \^ v }, (21a) 

(<S> p v l\H N + {H^T} + {H N S}\<S> V ) 

= E^(^Z\S2\^v)- (21b) 

Similarly, using the same definitions, the linearized equa- 
tion of triple excitation cluster operators is 

WW} + {H^S}\$ V ) = E«*(*$£\S 3 \Z V ).(22) 

The above equation of 5*3, except for the absence of H^, 
is very similar to the single and double excitation cluster 
equations. This key difference is due to structure of H^, 
which is either one- or two-body operator. The triples 
excitation operator S3 is, however, a three-body opera- 
tor. 

A more illustrative way to write the RCC equations is 
to identify the unique diagrams from the contractions and 
write the equivalent algebraic expressions. The linearized 
CC equation of singles, Eq. (21a), in terms of the cluster 
amplitudes is then 



(AE V - el)sl = 

bq 



ybp.q , -bp qr 

V qv L b ^ 2—1 q r bv 



E 

bcqr 



+ (23) 



IV. TRIPLES FROM LINEARISED RCC 

Extending the CCSD approximation in RCC to in- 
clude triple excitation is not a difficult proposition but 
entails enormous computational complications. In ad- 
dition, there is several orders of magnitude increase in 
the number of cluster amplitudes. The diagrammatic 
analysis, albeit easier and tractable, is cumbersome as 



Here, AE V — E% u — e v is the valence correlation energy 



of *„>,eg 



and xM = Vil — = — v l £ , is the 



e v - e p and v% = v% - WJ| = v l} ^ . 
antysymmetrized matrix element. Similarly, for compact 
notations the antysymmetrized closed-shell and valence 
cluster amplitudes are defined as ffj and ■ , respectively. 
Like in Si equations, we retain the T and S from the 
CCSD equations in S 2 equation as well. However, from 
triples cluster amplitudes, we consider only the valence 
triples S3. The Eq. (21b) in terms of cluster amplitudes 
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(A^ - <«K« = vg + E + «) - E «5!*E 

r c 

4- (v pC F q 4 7i c9 ;P r - 7 ,P c s r « _ ,;9 c f r P - 7 ; c P/ r « 

i / y \ u vr b cb "rb^vc u rb b vc u rv b bc u vr l cb 

rc 

- 11 CQ s r A 4 "i*" 8 4 Y^ 7J Cd / M 4 71 Cd f <t Prq 

v br s cv J T" 2^ u rs s i)fc + 2^ cd + 2^ \^udc 
rs cd cdr 

i s r P9 _ J>rq\ , ST^ (?pc srq _ pc srq\ 
' b vcd b vcd J ' Z^/ \ sr b vcb u sr b vbc ) 
src 

(24) 



P Q 
v <H- b 



Here, e^' = e„ 4 e& — e p — e q and ( ) indicates the 
combined permutations p f-> g and u -H- 6 of the previous 
terms within parenthesis. Interestingly, in this case terms 
with the combined permutations represent topologically 
distinct diagrams. For S3, the equation in terms of clus- 
ter amplitudes is 



(AE V 



vbc 



)s p 



be 



E 



E(„dr„pq ,,dr,pq n .pd.qr 
[ V bc S vd V vc Z db V vb T dc 



(„,qr „P S I „,pr „sq , pqj.sr\ 
[ V sc S vb + V sc S vb + V vs T bc J 



V. 



dr pqs 
bs h vdc 



rd pqs 
u sc b vbd 



„rd pqs 
J cs A vdb 



E' 

su 



,V1 ■ 

s u ' 



,,dr pqs 
u sc*vbd 



'vbc 



-E 

ds 
qd spr 
^bs b vdc 



q <-> r 
b o c 



- US 



_jpus 



V, 



de pqr 
be *vde- 



(25) 



Here, as defined earlier in the description of perturbed 
*^3' e «bc = e " + e b + £ c ^ e P — e 9 — In the present case, 
the combined permutations ( \ < ^ r c ) are just that, inter- 
change of the orbital lines and do not represent unique 
diagrams. Reason is, the two permutations q f-> r and 
b <-> c are between orbitals of the same kind virtual and 
core, respectively. Where as in S 2 equations one of the 
permutations is between core and valence, which have 
different topological representations. 



V. HFS CONSTANTS FROM RCC 

The hyperfine interactions Hm s are the coupling be- 
tween nuclear electromagnetic moments and electromag- 
netic fields of atomic electrons. The interaction energies 
from -ffhfs are the leading order corrections to the atomic 
and ionic energies obtained from H . In terms of the 
tensor operators, the hyperfine interaction Hamiltonian 

is mi es>] 



^h f8 = EE(- 1 ) 9 *"^) T - 9 ' 

i k,q 



(26) 



where vl{r) and T* are irreducible tensor operators of 
rank k in the electron and nuclear spaces respectively. 



For k = 1, following parity selection rules, the allowed 
interaction is the magnetic dipole. The explicit form of 
the associated tensor operators are 



ml _ 

1 q ~ /-*<?' 



-iy/2[a • Ci(f) 



(27a) 
(27b) 



where, Ci(f) is a rank one tensor operator in electron 
space and \x q is a component of fi, the nuclear magnetic 
moment operator. Interactions of higher rank multipoles 
are defined with similar form of tensor operators. How- 
ever, these are not discussed as in this work as we exam- 
ine the corrections to magnetic dipole hyperfine constants 
from the triples. From the expression in Eq. 
write the magnetic dipole HFS constant as 



2G 



we can 



A = 



(ttJ£i£ g (-lK(ri)Ti g l*t,) 



(28) 



where q = —1,0,1. The matrix element is calculated 
from the single particle reduced matrix element 



-.(nvKyW^WrivKv). (29) 



gi^N 

a = 

VJvUv + l)(2jv 

Here, gi (// = gil^N) is the gyromagnetic ratio, [in is the 
nuclear magneton and \n v K v ) is the valence single particle 
wave function. In a similar way, the HFS constants of 
higher order moments may be calculated. 
Using CC wave function from Eq. ^ 

(*«|ifhft|* ) - (*„|e Tt (l + 5) t ^hfiie T (H-S)|* ) J 

= (^l^hfe + 2S^H hfs 4 S^H his S\<I> v }.(30) 

Where, -ffhfs, = e Tt ^hfse T , is the dressed hyperfine inter- 
action and it is a non terminating series of closed-shell 
CC operator T. Further more, S^i/hfs = H hfs S is consid- 
ered while writing the equation. The higher order terms 
beyond second-order are, however, negligible and a trun- 
cated expression is considered. The approximation 



(31) 



accounts for all the important correlation effects and used 
in the present work. The normalization factor, denomi- 
nator in Eq. (28), is 



= (*„| (1 4 St) e Tt e T (1 4 SO I**). (32) 

Like in the dressed properties operator, e T ^e T is a non- 
terminating series. However, it is sufficient and accurate 
to consider up to the second order 

(*«|*»> w ($«| (l + S t S4T t T4S t T4T t S) |$ t ,).(33) 

The last two terms, although finite, are expected to be 
small as the contribution is of the form S|Ti and T^S2, 
respectively. For this reason, these two terms are not 
included in our calculations. 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



FIG. 5. Goldstone HFS diagrams from S2 fthfsS'3(vs), which 
contribute to A, vs y In the first diagram (a) all the two body 

vertices, Si, I/V12 and S2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



6 




(e) (f) (g) (h) 



FIG. 6. Goldstone HFS diagrams from S2 /ihf s 5 , 3(vs) , which 
contribute to Af va y In the first diagram (a) all the two body 

vertices, S%, l/ria and S2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



VI. HFS CONSTANTS FROM S 2 PERTURBED 
TRIPLES 



From the expression of HFS constant with RCC wave 
function in Eq. ( 30 1 , the lowest order triples contribu- 



tions are of the form 



A, 



hfe^S 



(34) 



The first term is, however, neglected in the present calcu- 
lations. The reason is, the T\ cluster amplitudes are small 
and have no significant contributions. For easy book 
keeping, contributions from the remaining two terms is 
bifurcated based on the nature of Hhfs matrix elements. 
Two of the possibilities, (a|/ih&|p) an d (w|/ihfe|p) are con- 
sidered. There are 52 Goldstone HFS diagrams asso- 
ciated with these two matrix elements and the S-x per- 
turbed S3. These are separated into groups and discussed 
in this section. The other forms, (a|/ihf s |a) and (p|fthfe|p)j 
enter through the structural radiation diagrams, which 
are negligibly small and are excluded from the present 
calculations. 



A. Contribution from S^H-yssSzt 



vs) 



Consider the triples of the form s^* 6 ^ 



defined in Eq. 

(12 1, for easy reference S^ vs ) define the general form of 

the triples in this group. The contraction /ihfs>53(vs) is di- 
agrammatically realized through four unique topologies. 
First, take the case where hhfs has a core-particle con- 
traction with S*2 in S^ vs ) and contribution to A is 

A 1 — \^ \^ \ ( „pq* ur pq sr _ pq* ir pq rs 
(vs) — 2^1 2-*i e P<l r va ' l b U sa b vb b va n b u sa b vb 
vab pqrs va b 

S va n b V sa S vb + S va rL b V sa S vb) i \ M ) 

where h T h denotes the matrix element (r|/ihf s |6). The 
many-body diagrams in Fig. [5^,-d are the representation 



of the above terms. For compact notation, introduce the 
antisymmetrised representation of the residual Coulomb 



matrix element, v p „i 



Antisym- 



metrised representation of the is defined in the same 
way. In a more compact form 



1 



jpqr °va IL b u sa°vb 



(36) 



vab pqrs 



It must be noted that, the antisymmetrised form is em- 
ployed for compact notations. Otherwise, all the calcu- 
lations are in non symmetrised representations and is a 
better choice with diagrammatic analysis. 

Second, the core and particle lines of ft-hfs contracts 
with the residual Coulomb and S2, respectively. Dia- 
grams arising from the contractions are shown in Fig. 
[5fe-h and contribution is 



(vs) — pqr \ b bv n a v sa b bv 



q r q*hP,, r i °P S 

b bv n a u sa b bv 



vab pqrs VOj b 

_J-q*ip qr sp , rq* , p qr ps\ 
b bv n a v sa b bv ' b bv n a v sa b bv)- 

In antisymmetrised representation 



A 2 



-1 



EE, 

vab pqrs va b 



„ r q*hPi, r 1 Z S P 

.pqr s bv ' L a v sa s bV 



(37) 



(38) 



The two cases discussed so far have double virtual orbital 
contraction of S\ with either residual Coulomb interac- 
tion or 5*2. As a result no unique diagrams arise from the 
anti-symmetrization of the S\ ■ 

Third, the core-virtual orbital lines of contracts 
with the residual Coulomb interaction. Eight unique dia- 
grams arise from the contractions and are shown in Fig|6] 
The contribution is 

A 3 — \^ \^ 1 \aPq*h b v qr *P s - q 9P*/ J b ,~,9 r sP ;s l 

(vs) — / , / , ^pqr l°va ll r u sb va °va "'r u sb va] ) 
vab pqrs va b 

(39) 
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(a) (b) (c) (d) 




(e) (0 (g) (h) 



FIG. 7. Goldstone HFS diagrams from S2 hhisSs( VB ), which 
contribute to A( vs -j. In the first diagram (a) all the two body 

vertices, Si, l/r-12 and S2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



FIG. 8. Goldstone HFS diagrams from S2 fthf s S , 3( cs ), which 
contribute to ^4( cs ). In the first diagram (a) all the two body 

vertices, Si, l/ri2 and 5*2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



here, as in previous expressions the antisymmetrised v 
and S2 are used for compact notations. The antisym- 
metrised expression of S% can be used to obtain the ex- 
pression 

Afvs) = E E ^Mr^vT^r^lb^va- (40) 
vab pqrs vab 

There is a prominent difference of the present case from 
the previous two, the exchange of S% gives topologically 
unique diagrams. 

Finaly, the core and virtual orbital line of /ihfs contract 
with the S2 and residual Coulomb interaction, respec- 
tively. Diagrams from the contractions are shown in Fig. 
[6] and in antisymmetrised notations, the contribution is 

A U = J212^%%W«- (4i) 

vab pqrs va ^ 

In this case too, there are eight unique diagrams. 



and are shown in Fig. [9] The contribution is 

A 2 \ \ ^ n pq* }J> ~cr ~pq l ao\ 

( cs ) _ / j / j ^pqr°va IL r v ab ; 'vc- \^°) 
vab pqrs va ^ 



Note, the exchange at the Si, like in -A( vs ) and A?s, does 
not generate topologically unique diagrams. 




(a) (b) (c) (d) 



FIG. 9. Goldstone HFS diagrams from /ihfsS^cs), which 
contribute to A^ CB y In the first diagram (a) all the two body 

vertices, Si, l/r\2 and 5*2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



B. Contribution from S^H^S^s) 



The triples of the <S3( CS ) type have two virtual lines 
above the S% vertex and no core line. This limits the 
number of allowed contractions between h^s and 5*2. So, 
there are only two unique topologies of the contraction 

^hfsS^cs). First, the core and virtual orbitals of h\^ s con- 
tract wit the residual Coulomb interaction and 5*2, re- 
spectively. Diagrams arising from the contractions are 
shown in Fig. [8] and contribution in antisymmetrised no- 
tation is 



A 1 



-1 



vab pqrs va ° 



c pqr J vb 



~pr* h a~cr~pq 



(42) 



And second, the h^s contracts with the orbital lines of 
residual Coulomb interactions. There are four diagrams 



C. Contribution from T%H Ms S 3 



There is a key topological difference between the 
r|_ffi 1 f s S f 3 and S , |_ffi 1 f s S , 3 diagrams. This arises from the 
number of lines above vertex of the cluster operators S2 
and Ti- The former has three, where as the later has four 
and more operators to contract. Consequently, fewer di- 
agrams arise from T^H^Sj, and these, like earlier, are 
identified based on the topology of contractions. Contri- 
butions from this term, like in S l Ji?hfs<S'3, is separable into 
7 1 |i?hts'S'3(vs) and T^H^S^^y Consider the first term, 
there are two groups of diagrams. In the first group, T 2 ^ 
contracts with the a pair of core and virtual lines with Si 



and v. Eight distinct diagrams, shown in Fig. 10 arise 



from this contraction and contribution is 



A 5 

(vs) 



EE, 

vab pqrs va ° 



—t qr *h v v pr s sq (W) 

c pqr b ab n p u sb°va- V**J 



Here, we have given the antisymmetrised expression. The 
individual terms may be written in explicit forms line in 
Eq. (39). The second group of diagrams arise from the 




FIG. 10. Goldstone HFS diagrams from S2 frhfsiS3(vs), which 
contribute to -A( vs \. In the first diagram (a) all the two body 

vertices, Si, l/ri2 and S2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



contraction of T2 with two virtual and one core orbital 
lines of v, and one core orbital line of 6*2. Four diagrams 
arise from this term and these are given in Fig. [TT] The 
contribution in antisymmetrised form is 



A 6 



EE e . 

vab pqrs 



_fsr* h v sr~pq (az\ 
pqr b ab "'p' J qb c 'va- V*"; 
vab 



Note that T2 and S2 are antisymmetrised in the above 
equation. Where as, v and S2 are antisymmetrised in the 
previous groups consisting of four diagrams. These two 
antisymmetrizations are equivalent and give the same set 
of diagrams. The term ^4( vs ) completes the possible forms 
of HFS diagrams arising from the <S3(vs) type of valence 
triples. Collecting all the terms, the net contribution is 




(d) 



FIG. 11. Goldstone HFS diagrams from S2 frhfsSWvs), which 
contribute to A^ VB y In the first diagram (a) all the two body 

vertices, S|, I/V12 and 5*2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



The second group has two diagrams and these are shown 
in Fig. [l2fe-f. The contribution is 



(cs) — f P<ir l ab n r v ab b 



ZP1 
dv " 



(48) 



vabc pqr va ^> 



Collecting all the groups, the net contribution from the 
SWcs) type of triples is 



.4 



(cs) 



4 

E A 



(cs)- 



(49) 



i=X 



Totally there are 18 Goldstone diagrams in Ar cs y Col- 
lecting all the diagrams from the S2 perturbed triples, 
we define 



A(s) — A( vs ) + A( cs y 



(50) 



There are in 54 many-body diagrams and these are sep- 
arable into ten groups. This completes, excluding the 
structural radiation diagrams, the diagrammatic analy- 
sis of the S2 perturbed triple correction to the magnetic 
hyperfine constant. 




(b) 




A 



(vs; 



= E^ 



(vs) 



(46) 



(e) 




To summarize, Ar vs \ constitute 36 many-body Goldstone 
diagrams grouped into six groups. Each group is defined 
based on the contraction topology and form of ft.hf s - 

From r^i/ h fs5 , 3(cs) there are two groups of diagrams. FIG. 12. Goldstone HFS diagrams from S 2 t /ihfeS , 3 ( cs ), which 

The first group has four diagrams and these are shown contribute to Af cs) . In the first diagram (a) all the two body 

in Fig. [T2fc,-d. The contribution is vertices, Si , 1 /r 12 and S 2 are of direct type. Remaining dia- 
grams are combinations of exchange at different vertices. 



A 3 



EE 

vabc pqr 



1 rqr*,v cr~pq 
.pqr L ab IL p u ab°vc- 



(47) 
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VII. HFS CONSTANTS FROM T 2 PERTURBED 
TRIPLES 



The contributions from the T 2 perturbed triples, like 
in S2 perturbed triples, is separated into two categories: 
perturbation to the core orbital and virtual orbital. Con- 
tribution from these are defined as -A( ct ) and ^4(vt) ■ Simi- 
lar to Ar a \ , diagrams from each of these are classified into 
groups. Diagrams from each of the groups with direct at 
all the two body vertices are given in Fig. [l3j The dia- 
grams of the -A(ct) are shown in Fig. 13'a-d) and (g-h). 
The expression is 



Act) - X! X! ~^pW 

vabc pqr va ^ 



,qp* ub-.crfqp _ „qr*ub~cplqr 
'av "'r u vb b ac °av "'p u av L cb 



+s„ h n T v„ h x nh + s vb n q v va i cb -t i ab n p v vb i ac 



vb r u vb cb 

+t qr *h v v pc t qr 

^ b ab n p u va L cb 



(51) 



One immediate observation is, the structure and number 
of the terms in the above expression are similar to At vs -\ . 
Key transformations are interchange of conversion of S2 
and v operators to v and T2, respectively. The expression 
of the second category is 



TABLE I. Basis set parameters 


a and ft used 


in the calcula- 


tions. 










Atom 


Orbital 


a 


P 


Basis function 


Rb 


s 


0.00521 


2.9500 


33 




P 


0.00655 


2.9950 


30 




d 


0.00654 


2.9720 


28 


Sr+ 


s 


0.00825 


2.9000 


35 




p 


0.00715 


2.9450 


32 




(1 


0.00730 


2.9100 


30 



VIII. RESULTS AND DISCUSSIONS 

A. Single particle states 

The first step of our calculations, like in any atomic 
many-body calculations, is to solve the single particle 
eigenvalue equations with Dirac-Hartree-Fock potential. 
For this, we consider the nuclear potential V^(r) arising 
from the finite size Fermi density distribution 



3 (r) 



Po 



1 + e^-c)/ ' 



(53) 



Avt) - X! e pi r 

vab pqrs va ^ 



vb "'q v vs L ab 



s P1*h b yPlf Sr ; 
°va " , r v vs l, ab 



+i q ; b *h;z p v q s t% - tz*Knrx r b ) ■ (52) 

Here, the terms are similar to A^ and same transfor- 



mations discussed in ^4( c t) apply- 




FIG. 13. HFS diagrams from the T2 perturbed 5*3. Here, 
only the diagrams with direct interaction at all the interaction 
vertices are given. 



here, a = £4 In 3. The parameter c is the half-charge 
radius, that is p n uc(c) = po/2 and t is the skin thickness. 
At the single particle level, the spin orbitals are of the 
form 



1 



r \ iQ nK (r)x-« 



>(?) 
>(r) 



(54) 



where P raK (r) and Q nK ,{r) are the large and small compo- 
nent radial wave functions, k is the relativistic total an- 
gular momentum quantum number and Xumi?) are the 
spin or spherical harmonics. One representation of the 
radial components is to define these as linear combination 
of Gaussian like functions and are referred to as Gaus- 
sian type orbitals (GTOs). Then, the large and small 
components [301131] are 



^(r) = DC^(r), 
p 



(55) 



The index p varies over the number of the basis functions. 
For large component we choose 



(r) 



(56) 



here n K is an integer. Similarly, the small component are 
derived from the large components using kinetic balance 
condition. The exponents in the above expression follow 
the general relation 



a /3 : 



(57) 
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The parameters ao and f3 are optimized for each of the 
ions to provide good description of the properties. In 
our case the optimization is to reproduce the numerical 
result of the total and orbital energies. The optimized 
parameters used in the calculations arc listed in Tablc|T] 



TABLE II. Excitation energies calculated using RCC, com- 
pared with other theoretical results and experimental data. 
All values are in atomic units. 



Atom 


State 


This work 


Other works 


Exp Ref[55], 


S5 Rb 


5Sl/2 


0.0 


0.0 


0.0 




5pi/2 


0.05759 


0.05718 a 


0.05731 




5P3/2 


0.05872 


0.05826 a 


0.05840 




4d 3 /2 


0.08836 


0.08822 a 


0.08819 




4d 5 /2 


0.08836 


0.08820 a 


0.08819 


87 Sr+ 


5S1/2 


0.0 


0.0 


0.0 




5pi/2 


0.10841 


0.11001 b 


0.10805 




5P3/2 


0.11221 


0.11376 b 


0.11171 




4d 3 /2 


0.06611 


0.06560 b 


0.06632 




4d 5 /2 


0.06724 


0.06707 b 


0.06760 



Reference [34] Reference [35] 



For Rb and Sr+ we use V N_1 and V N ~ 2 orbitals, re- 
spectively. These are the single particle eigenfunctions of 
the Rb + and Sr 2+ ions, respectively. The single particle 
basis sets have few bound states and rest are continuum. 
We optimize the basis such that: single particle energies 
of the core and valence orbitals are in good agreement 
with the numerical results. For this we use GRASP92 
[3"2"] to generate the numerical results. It is to be noted 
that, the basis parameters in Table|l] are different from 
one give in our earlier work [24]. Between the two, the 
present is better optimized and of higher quality. The 
optimization is nontrivial as there are several parameters 
and single particle equations are solved self-consistently. 

From Eq. ( 54 1 the reduced matrix element of the mag- 
netic hyperfine operator between two spin orbitals , v' 
and v, is 



(t/|t 1 |v) = ~(k v + K t /)(-K l ,'||C 1 ||K„) 
CLT 
r 



X 



A detailed derivation is given in Ref. 



B. Cluster amplitudes and normalization 



As described in our earlier works [23, 24, 36 , the RCC 
equations are solved iteratively using the Jacobi method. 
To improve convergence we employ direct inversion in the 
iterated subspacc (DIIS) [37]. The cluster amplitudes 
are solved for each Hilbert space manifold of the total 
Fock space. At each step the Hilbert space is augmented 
with one electron. In short, the T equations are solved 



first and these are used to generate the open shell cluster 
amplitudes 5*. 

One important point is, the cluster equations are in 
terms of the reduced matrix elements. So the solutions 
are independent of magnetic quantum numbers and ap- 
propriate phase factors are required to define the cluster 
amplitudes in the cojugate manifold and these are 



a 6 



(-i) jW -e 



(59) 
(60) 



These relations apply in any calculation which involve 
and . The coupled-cluster wave function is normalized 
and the normalization factor is 

N = (¥„|¥ ) = ($„|e Tt (1 + 5+) (1 + S) e T |$„). (61) 

Here, e Tt e T is a non-terminating operator. For the 
present we consider the approximation 

N fa ($t>|l + SflS + T^T\$ V ). (62) 

The higher order terms (T 1 t ) a (T 2 t ) /3 (T 1 ) 7 (T 2 )' 5 , such that 
a + fi >1, 7 + 5>l and a + (i — 7 — 5 — 0, are neglected. 
We also neglect the mixed operator term (T^S^ + c - c 
and higher orders. 



C. Excitation energies 

To determine the quality of the basis set and parame- 
ters, we compute the attachment energies of the ground 
state (S1/2) and the first excited P1/2, P3/2, D3/2 an d 
-D5/2 states are calculated. Then the ionization potential 
(IP), the energy required to remove the valence electron, 
is the negative of the attachement energy — E att . To 
calculate the excitation energy (EE) of the state \^ v ), 
consider E att and E* tt as the attachment energies of 
the ground state and excited state. Then difference 



of IPs. 



Eg U is the EE, it can as well be defined in terms 



D. HFS constants 

To compute the hyperfine constants from the CCSD 
wave functions, we use Eq.(30|. The results are listed 



in Table III for comparison the results of other theoret- 
ical calculations and experimental data are also given. 
As defined in Eq.(30l, the coupled-cluster expression of 



the hyperfine structure constants is separated into three 
groups. The dominant contribution from the first term 
Hhfs, up to first order in and T, is 



Hhfs ~ -ffhfs + 2i7hf s Ti + T x Hhfs (Ti + 2T 2 ) + T 2 ifhfsT^. 

(63) 

Here, the first term is the Dirac-Fock (DF), which has the 
largest contribution. The factor two in the second and 
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TABLE III. Magnetic dipole 


hyperfine structure constants for 


85 Rb and 87 Sr+. The values given 


are in the unit of MHz. 


Atom 


State 


This work 


v./ LUcI VVUi K.Q 








CCSD 


CCSD(T) 






85 Rb 


5S1/2 


1030.94 


1030.60 


1011 1 a 


1011 Q1 (9\ b 




5pi/ 2 


120.69 


120.67 




1 90 79('9^ 1 \ b 




5P3/2 


24.48 


24.63 


24. 5 a 


94 QQfl ~l c 




4053/2 


7.85 


7.80 




7 ^f^WW 7 ^9Qf*3^ e 




4rf 5 /2 


-4.78 


-4.77 


_ 


-5.2(3) d 


87 Sr+ 


5S1/2 


-1014.20 


-1013.86 


-10003.18 f ,-1000 g 


1000.5(1.0) j 




5P1/2 


-178.73 


-178.67 


-178.40 f ,-177 s ,-175 h 






5P3/2 


-35.28 


-35.48 


-35.11 f ,-35.3 g ,-30 h 


-36.0 j 




4d 3 / 2 


-46.30 


-46.21 


-47.36 f ,-46.7 g 






4d 5/2 


1.71 


1.71 


2.51 f ,l.l g ,2.156 i 


2.17 k 


a Reference [55], b Reference [3"9] 
h Reference [45] . 1 Reference |46|. 


c Reference [4*0] 
j Reference [4"7]. 


d Reference [H], 
k Reference [48], 


e Reference [12], f Reference gS], g 


Reference [4"4"]. 



TABLE IV. The HFS contributions from different terms in RCC. All the values are in MHz. 



Atom 


State 








RCC terms 














DF 


iihfe-DF 


S^H his 


SlHhisSl 


S 1 HhlaSl 


S 2 HhlaS2 


Other 


Norm 










+C.C 


+C.C. 










85 Rb 


5«l/2 


646.003 


-12.115 


375.719 


19.822 


23.564 


20.498 


-22.231 


1.020 




5P1/2 


69.883 


-0.762 


47.411 


2.762 


2.751 


1.630 


-1.760 


1.010 




5P3/2 


12.417 


-0.054 


10.919 


0.804 


0.450 


0.479 


-0.312 


1.009 




4d 3 /2 


3.116 


0.127 


3.153 


0.119 


1.014 


0.685 


-0.068 


1.039 




4d5/2 


1.328 


0.087 


-4.694 


-2.550 


0.420 


0.472 


-.028 


1.038 


87 Sr+ 


5*1/2 


-735.629 


8.044 


-284.856 


-9.987 


-8.594 


-15.624 


16.907 


1.015 




5P1/2 


-122.123 


1.369 


-56.502 


-2.155 


-1.862 


-1.566 


2.217 


1.010 




5P3/2 


-21.449 


0.191 


-13.023 


-0.582 


-0.307 


-0.869 


0.412 


1.010 




4d 3 /2 


-31.401 


-0.803 


-12.400 


-0.594 


-0.428 


-2.020 


0.283 


1.023 




4d 5/2 


-13.091 


-0.428 


16.306 


1.051 


-0.177 


-2.018 


0.116 


1.023 



fourth terms accounts for the complex conjugate terms. 
The third term, second order in T\, has one diagram and 
negligibly small contribution. The diagrams arising from 
the last term are topologically are the structural radi- 
ation diagrams and have neglible contributions. Topo- 
logically, these are insertion of H^ s to the normaliza- 
tion diagrams and contribution from these are labelled 
as Hm s — DF. Detailed diagrammatic analysis are given 
in our previous work [24]. The last two terms in Eq.(30) 
are approximated as 



&H Ms w 25+ (H hfs e q 
S'HbfsS w S 1 Hms (Si 



2S 2 



5|i/hf s 52. 



(64) 
(65) 



Like in i/hfs, the factor of two is to account for the 
complex conjugate terms. Based on this grouping, the 
contributions are listed in Table |TV] In the following we 
present a detailed comparison of our magnetic hyperfine 
constants results with the earlier ones. As discussed later, 
some of our results are the best match with experimental 
data. 



E. HFS constants contribution from triples 



The HFS constants after including the perturbed 
triples are listed in the Table. 



Ill 



There is negligible 
contribution for Sin and Pin states. These are 0.03% 
and 0.02% for Rb and, 0.03% and 0.03% for Sr+. How- 
ever, for P3/2J -D3/2) & n d -D5/2 contributions from triples 
are not small and could be important in high precision 
atomic theory calculation. And these are 0.6%, 0.6% 
and 0.2% for Rb, and 0.6%, 0.2% and 0.5% for Sr+. The 
observed pattern of perturbed triples contribution is dif- 
ferent from 87 Rb reported in Ref. [51]. This could be on 
account of two factors: difference in the nature of the sin- 
gle particle basis functions, and isotope specific effects. 
The later may not be the dominant cause as the electron 
wavefunctions have little variation for isotopes of small 
mass differences. For Sr + on the other hand, there are 
no previous theoretical work on the effects of triple ex- 
citations. However, our results exhibits trends similar 
to previous work on Ca + [25] , which reported the con- 
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tributions from triples as 0.002%, 0.08%, 0.10%, 0.11% 
and 0.29% for 5i/ 2 , Pi/ 2 , -P3/2, -D3/2, and D 5/2 states, 
respectively. 



TABLE V. The term wise contribution of HFS constants from 
triples for Rb. All values listed are in atomic units. 



RCC term 



A 2 

A (y 

A 3 
A (v 

A° 



5fi/2 5pi/ 2 5p 3/2 



4d 



3/2 



4rf, 



( vs 



T 



(vs 



A \ 
4 

A? 



-4, 



r 

(ct 
5 

1 (ct 
(ct 



A 2 



A 



ill 



(vt 



0.083 
0.027 
0.200 
0.064 
-0.236 
-0.135 



0.013 
0.000 
0.003 
-0.002 
-0.017 
-0.014 



0.013 
0.014 
0.024 
0.012 
-0.010 
-0.008 



-0.003 
0.004 
0.020 
0.000 
-0.002 
-0.005 



-0.054 
0.002 
0.000 
-0.028 
-0.002 
-0.004 



0.003 

-0.141 
0.012 
0.216 

-0.037 



-0.017 

-0.007 
0.007 
0.018 

-0.002 



0.045 

0.024 
0.000 
0.010 
-0.001 



0.014 

-0.018 
0.001 
0.007 

-0.001 



-0.086 

0.058 
0.015 
0.005 
-0.001 



0.050 

-0.219 
-0.007 
-0.005 
-0.007 
-0.149 
1.966 
1.579 

-0.084 
0.008 

-1.782 
0.054 



0.016 

-0.006 
-0.002 
0.000 
-0.002 
-0.015 
0.195 
0.170 

-0.006 
-0.001 
-0.178 
0.006 



0.033 

-0.005 
0.000 
-0.002 
-0.002 
-0.007 
0.113 
0.097 

0.002 
-0.002 
-0.099 

0.004 



-0.011 

0.000 
-0.005 

0.001 
-0.002 
-0.011 

0.035 

0.018 

-0.014 
-0.003 
-0.028 
0.004 



0.077 

0.000 
-0.001 
-0.003 
-0.007 
-0.007 
0.030 
0.012 

0.007 
0.014 
-0.023 
0.003 



-1.804 



-0.179 



-0.095 



-0.041 



0.001 



Rb 



In the Table. |Vj we have listed the individual contri- 
butions from the different groups of triples HFS diagrams 
for Rb. In all the cases we notice large cancellations. For 
S1/2 and P 1 / 2 states, the leading order (LO) and next 
to leading order (NLO) contribution arise from A® t and 
A^ t . Each of these terms contribute pa 0.2%. However, 
the two are of opposite signs and cancel each other. The 
other dominant contributing terms are A^ s and A^ s with 
contrbutions pa 0.02%. The two contributions are op- 
posite in sign and like in LO and NLO, there arc large 
cancellations. 

For the state P3/2, A® t and A^ t are again the LO and 
NLO terms with contributions of pa 0.5% and pa 0.4%, 
respectively. The two contributions are of opposite sign 
and nearly cancel. Other dominant contributions arise 
from A^s and A 1 ., each of the contributions are pa 0.1%. 
Unlike the cases considered and discussed so far, the two 



TABLE VI. The term wise contribution of HFS constants 
from triples for Sr + . All values listed are in atomic units. 



RCC term 



~A 1 ~ 

A lvs 

A 2 

7 s 

A (va 

A% B 

A 5 
A( V s 

A 1 

A 2 

A A 



5fl/2 5pi/ 2 5p 3 / 2 



Ad, 



3/2 



4d r , 



11. 



A( C s 

Al t 
AL 

a) 

A A 
A^ 



A 



(ct 



Al t 
A? 

A (vt 

aU 



A 



(vt 



-0.066 
-0.020 
-0.129 
-0.045 
0.212 
0.109 
0.061 

0.092 
-0.010 
-0.157 

0.028 



-0.016 
0.001 
0.001 
0.000 
0.027 
0.018 
0.031 

0.008 
-0.004 
-0.022 

0.003 



-0.015 
-0.020 
-0.039 
-0.015 
0.016 
0.011 
-0.062 

-0.028 
0.000 

-0.012 
0.003 



-0.011 
-0.005 
-0.046 
-0.001 
0.012 
0.026 
-0.025 

0.036 
-0.010 
-0.032 

0.006 



0.098 
0.001 
0.024 
0.070 
0.010 
0.020 
0.223 

-0.098 
-0.055 
-0.024 
0.004 



-0.047 

0.158 
0.006 
0.003 
0.005 
0.107 
-1.935 



-0.015 

0.007 
0.004 
0.000 
0.002 
0.017 
-0.302 



-0.037 

0.008 
0.002 
0.004 
0.001 
0.008 
-0.172 



0.000 -0.173 



0.000 
0.016 

-0.001 
0.005 
0.092 

-0.327 



0.000 
0.002 
0.006 
0.024 
0.059 
-0.266 



-1.656 

0.069 
-0.003 

1.786 
-0.041 

1.811 



-0.272 

0.010 
0.001 
0.281 
-0.008 
0.284 



-0.149 

-0.003 
0.002 
0.155 

-0.004 
0.150 



-0.215 

0.038 
0.010 
0.268 
-0.033 
0.283 



-0.175 

-0.030 
-0.032 

0.209 
-0.025 

0.122 



P3/2, A® t and A^ t are the LO and NLO terms. Contribu- 
tions from these terms are ps 0.4%, but opposite in sign. 
The other dominant contributions are from A^ s and Aj s 
and these are « 0.3% and pa 0.2%, respectively. Like in 
LO and NLO these are of opposite sign and nearly cancel. 

The state D 5 / 2 shows a different pattern of contribu- 
tions. Unlike the states discussed so far, the dominant 
LO term is A 1 , and the contribution is about 1.2%. The 
NLO term is A* a and has a contribution of psI.1%. and 
is opposite to the LO term. The terms A® t and A^ t , 
which are the LO and NLO of <Si/ 2 , P1/2: P3/2 and -D3/2, 
are the third and fourth dominant terms. Contributions 
from these terms are pa 0.6% and pa 0.5%, respectively. 



2. Sr+ 



The contributions from different RCC terms for Sr + 



are as listed in the Table. VI For the states S1/2, P1/2 



are of same phase. For -D3/2 too like in S l 1 / 2 , P1/2 an d 



/2- 



and -P3/2, the LO and NLO have the same pattern as Rb 
is observed. The only difference is, signs are reversed. 
For example, in the case of S\/2 of Rb, A^ s and A® t are 
positive in sign, and A^ s and A^ t have negative sign. 
The pattern is opposite in the case of Sr + . The reason is 
trivial and is on account of the opposite sign of fi, which 
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TABLE VII. Excitation energies calculated using RCC for 
different sets of basis functions. All values are in atomic units. 

Atom State Set 1 Set 2 Set 3 Set 4 Set 5 

"S3 - 



3 Rb 



7 Sr+ 



5Si/ 2 
5pi/2 
5P3/2 

Ad 
4d fi 



3/2 
/2 



5Si/ 2 
4^3/2 

4d 5 /2 

5P1/2 
5^3/2 



0.0 0.0 0.0 0.0 0.0 

0.05756 0.05754 0.05758 0.05759 0.05759 

0.05869 0.05867 0.05871 0.05872 0.05871 

0.08854 0.08854 0.08834 0.08836 0.08836 

0.08854 0.08853 0.08834 0.08836 0.08836 

0.0 0.0 0.0 0.0 0.0 

0.06772 0.06767 0.06611 0.06611 0.06607 

0.06882 0.06877 0.06724 0.06724 0.06720 

0.10836 0.10837 0.10840 0.10841 0.10841 

0.11217 0.11217 0.11220 0.11221 0.11221 



TABLE VIII. The HFS constants calculated using RCC for 
different sets of basis functions. All values listed are in units 
of MHz. 



Atom State 


Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


S5 Rb 5s 1/2 


1029.19 


1027.78 


1030.70 


1031.00 


1030.94 


5^1/2 


120.66 


120.59 


120.66 


120.71 


120.69 


5P3/2 


24.47 


24.45 


24.48 


24.48 


24.48 


4d 3 /2 


7.54 


7.52 


7.84 


7.84 


7.85 


4d 5 /2 


-4.78 


-4.74 


-4.79 


-4.78 


-4.78 



7 Sr+ 5si/ 2 
5pi/ 2 

5P3/2 

Ad 



3/2 
4d 5 /2 



-1012.82 -1012.97 -1014.09 -1014.20 -1014.20 
-178.61 -178.63 -178.67 -178.75 -178.73 
-35.25 -35.25 -35.28 -35.28 -35.28 
-45.90 -45.92 -46.21 -46.31 -46.30 
2.085 2.07 1.73 1.72 1.71 



for Rb and Sr + are positive and negative, respectively. 

For -D3/2 state, the LO and NLO contributions are 
large in comparison to the case of Rb. The LO and NLO 
terms are A® t and A^ t , respectively, and contributions are 
« 0.7% and « 0.6%. The other dominant contributions 
are from A^ s and A^ s , and these are smaller than Rb. 

The state -D 5 / 2 show the largest LO and NLO contri- 
butions among all the states. These are 15.5% and 12.2% 
from the terms A® t and A^ t , respectively. However, the 
total contribution is 3.3% as they are of opposite signs. 
The other dominant contribution arise from A^ s and Al s . 
Contributions from each of these terms are 5.7% but are 
opposite in sign. 



F. Uncertainty estimates 

Atomic properties calculated from RCC theory, in gen- 
eral, have three important sources of uncertanties. These 
are: omission of higher-Z orbital basis states, truncation 
of the dressed HFS operator iJhfs and truncation of the 
coupled cluster operator T. The error arising from the 
third source-is truncation of the CC operator-is almost 
mitigated with the inclusion of perturbative triples. So, 
effectively, HFS results presented in the Table. 



Ill have 



uncertainties from the first two sources. In the following 
we analyze and estimate the upper bound on the uncer- 
tanties arising from each of these sources. 

The results presented in Table. |III| are the converged 
results with the orbitals up to h symmetry. To define 
the converged basis set, we do a series of calculations 
where we start with a minimal basis size of 112 orbitals 
consisting of (l-14)s, {2-U)p, (3-15)d, (4-16)/, and (5- 
14)g. Where we have used the non-relativistic notations 
for compact representations. The basis set is increased 
by adding two orbitals in each symmetry in the succes- 
sive sets of calculations till the change in the excitation 
energies and HFS constants are below 10~ 4 . The values 
from different orbital basis sets are given in Table. |VII| 
and |VIIl| The total number of orbitals in the converged 
results is 177, and symmetry wise it is (l-19)s, (2-19)p, 
(3-20)d, (4-21)/, (5-19)3 and (6-15)/i. The single particle 



energies considered are « 10 5 Hartrees for s, p, d and / 
orbitals and, sa 10 3 Hartree for g and h. 

To estimate the uncertanties from excluding higher 
symmetries, we include orbitals of i symmetry and com- 
pute the HFS constants. For S1/2, the largest contribu- 
tion is in the case of Rb and is about 0.07%. However, 
for the states P1/2 and P 3 / 2 it is Sr + which has large 
contribution from the i symmetry. These are 0.03% and 
0.04% respectively for the states P1/2 and P 3 /2. Unlike 
the S1/2, P1/2 and P 3 / 2 states, the triples contributions 
for .D3/2 and -D5/2 are in general large. These are 1.2% 
for -D3/2 in the case of Rb and 7.0% for D5/2 in the case 
of Sr + . In the present implementation of RCC theory, 
incorporating j and higher symmetry basis states ren- 
ders the basis set too large for computations. However, 
a leading-order analysis is possible with MBPT and we 
find the contribution from j symmetry is negligible. To 
estimate the uncertanties arising from the truncation of 
dressed properties operator, we resort to our previous 
work |24j . Where we proposed an iterative scheme to 
account for the higher order terms in the dressed prop- 
erties operator Phfs- Using this scheme, we computed 
HFS constants contributions from the RCC terms which 
have loe of zero and one. These are the categories of di- 
agrams which contribute most to the atomic properties 
using RCC. Contributions from the loe zero are 0.009%, 
0.008%, 0.02%, 0.03% and 1.5%, respectively for S 1/2 , 
P\/2i P3/2: -D3/2 and -D5/2 states of Sr + . From the loe 
one, however, these are large 0.09%, 0.24%, 0.35%, 0.22% 
and 12.40%. The loe two and higher are not considered 
as it involve higher-order terms in T which will naturally 
have smaller contributions. 

In addition to the sources of errors discussed so far, the 
other sources of errors are the QED corrections. How- 
ever, this cannot be considered as the source of error of 
the many-body method. It is rather the form of interac- 
tion considered in the calculations. 

To put an upper bound on the uncertainty in the HFS 
results, we select and add the largest change from each of 
the sources. Which turn out approximately to be 0.2%, 
0.3%, 0.4% and 1.5% for the states S 1/2 , P1/2, P3/2 and 
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D3/2, respectively. For the state D 5 / 2 , since there is large 
cancellations, a comprehensive analysis is necessary to 
arrive at a meaningful uncertainty estimate. 



IX. CONCLUSIONS 

We derive and propose a general tensor representa- 
tion of the triple cluster operator is symmetric and with 
proper phase the core and virtual orbital indices may be 
permuted. The permuation properties are derived from 
the rules of angular momentum diagrams. Although, we 
have discussed the triple cluster operator in the context 
of valence triples S3, the same definition applies to core 
triples T 3 . This may be explicitly demonstrated with a 
minor topological transformation of the cluster operator. 
Based on the generalized tensor representation of 53, we 
derive the linearized RCC equation for CCSDT approxi- 
mation. 

We have analysed the contributions from the triple ex- 
citation cluster operators to the HFS constants in detail. 
For the analysis we identify two groups of diagrams de- 
pending on the nature of contractions to generate the 
triple cluster operator. In most of the cases the LO and 
NLO are identified as A% and Af rt . These arise from the 
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